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TECHNICAL MEMORANDUM X- 64594 


APPLICATION OF REGRESSION ANALYSIS TECHNIQUES 
TO REFRACTORY COATING MEASUREMENT 
EXPERIMENTS 


SUMMARY 


This report describes a procedure for conducting a statistical analysis of data obtained 
from two nondestructive techniques/instruments used to measure the thickness of protective 
refractory coatings on Columbium alloy. A regression analysis is performed on instrument 
output data to determine the significance of linear, quadratic, and cubic models. It is con- 
cluded that the linear regressions are the best choice of models for the particular coating, 
instrument, and base metals used. 


INTRODUCTION 


This report is the result of a request for computational support and analysis from 
the Materials Analysis Section of the Quality and Reliability Assurance Laboratory, Marshall 
Space Flight Center. A procedure is presented for conducting a statistical analysis of data 
obtained from two nondestructive techniques/instruments used to measure the thickness of 
protective refractory coatings on Columbium alloy. The two instruments used to obtain the 
test data are the Dermitron and the Betascope (a SPACO internal note gives a complete 
discussion of the test program and its relation to the Space Shuttle program 1 ). The evalua- 
tion and analysis of data obtained from these instruments are required to evaluate accepted 
and reliable thickness measuring principles. Implementation of the various statistical methods 
required in the evaluation and analysis effort is accomplished on the UNIVAC 1 108 com- 
puter. A complete program listing and input data setup examples are given in the appendix. 

Regression analyses were performed on the output data from both instruments. 
Linear, quadratic, and cubic models were assumed and specific hypotheses tested concerning 
the significance of the various models. Confidence limits were attached to the models to 
establish a range of accuracy at a given confidence level. It was found that the linear 
regressions are the best choice of models for the particular coating, instrument, and base 
metals used. The accepted models allow instrument response to be related to actual coating 
thickness. 


1 . Charles Wages and Marshall Parks: Evaluation of an Eddy Current and a Beta Backscatter 
Instrument for Measuring the Thickness of Refractory Coatings. Internal note, SPACO, 
Inc., Huntsville, Ala., Jan. 14, 1971. 



MATHEMATICAL DEVELOPMENT FOR MULTIPLE 
REGRESSION ANALYSIS 


General 

We assume that the observed response variable (dependent variable) Y- 3 is to be 
estimated by the model: 


Yf = b 0 + b 1 Z li + b 2 Z 2i + ... + b p Z pi 


pi 


( 1 ) 


where i = 1, 2, . . . , n. If the input data are centered about the mean, this model becomes: 


Y i * b 0 + b l z li + b 2 z 2i + 


+ b p z pi 


( 2 ) 


where 


z ii - z ir z i 


z 2i " Z 2i" Z 2 


z pi Z pi ” Z p 


(3) 


bo = b 0 + bj Z! + b 2 Z 2 + • • . + bp Zp 

The sample estimates b|> , b t , b 2 , . . . , b p are obtained by minimizing the weighted sum of 

squares of deviations between the observed and predicted response values [ 1 ] . It can be 
verified from the first normal equation in the least-squares formulation that bo — Y. This 
leads to the following equation: 


Yf = Y + b 1 z li + b 2 z 2i + ... + bpZ pi 


(4) 


2 



Consideration of the remaining normal equations leads to the following equation for the 
regression coefficients: 


B = s p S y (5) 

pX 1 pX p pX 1 


where 
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and the S’s are given by: 



% = 2 < z ii - Z I) czji - z j) 


S IY 



(Zji-ZjHYf-Y) 


I = 1 , 2, . . . , p for 
each J = 1,2, ... ,p 

I = 1, 2, . . . , p 


(9) 


The standard deviation of the dependent variable is given by: 


cry = 





n - d 


1/2 


( 10 ) 


where d = p + 1 . For the regression coefficients we have: 


a R “ °Y C (11) 

(p+l)Xl (p+l)Xl 

where 



a R 


(P+l)xi 


ff b 2 


a 



( 12 ) 
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(13) 


C 

(p+ 1 )X 1 


1/Vn 

V^n 


VC22 


Y^-pp 


The C’s are elements in the inverse of the S matrix as given by equation (7). 

Standard Error of Y c 

Consider equation (4). If Y, b, , b 2 , . . . , b p are subject to error and there is no 
correlation between Y, bj , b* , . . . , bp, the following error formula for the variance of 
Y c can be written: 



Confidence Bands 

The 95-percent confidence bands for the regression curve can now be established 
[2] by use of the Students T table. Let t T denote the appropriate table value for N 

degrees of freedom where : 

N = n-(p+ 1) ( 15 ) 

c 

Then the upper and lower bands for Yj are determined from: 


y[(U) = Y® + t T aye 
Y?(L) = Yf-t T a y c 


(16) 
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The resulting curves about the regression curve are the loci of the 95-percent confidence 
bands. These bands can be interpreted as follows. If repeated observations of y£ are 
taken and at the same fixed values of Zj^, Z2^, ■ . ., Zpj, as were used to determine the 
fitted regression equation, then of all the 95-percent confidence intervals constructed for 
the mean value of y£, 95 percent of these intervals will contain the true mean value of 
y£; or, there is a 0.95 probability that the following statement is correct: 


“The true mean value of y£ at Zj^, Z 2 jj, . . . , Zp^ lies in the interval 

[Yu •• <T »y>] •" 


Significance of the Estimated Regression Equation 

We consider the following analysis of variance table for the regression equation as 
given by equation (4). 


TABLE 1. ANALYSIS OF VARIANCE 


df 

Type 

Variation 

SS 

MS 

F Value 

n- 1 

Total 

S YY = t (Yf-Y) 2 

£l 



n - d 

Residual 

n 

S(RES) = ^ (Y? - Yf) Z 
f1 

M(RES) = SffiESl 
n - d 


d- 1 

Regression 

n 

S(REG) = ^ (Y? - Y) 2 

M(REG) = S(REG) 
d - 1 

M(REG) 

M(RES) 

df = degrees of freedom 

SS = sum of squares 

MS = mean square = SS/df 

n = number of observations 

p = number of independent variables 

d = p + 1 
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The total sum of squares can be written as: 


S Y y = S(RES) + S(REG) 


( 17 ) 


If S(RES) = 0, the actual observations of the dependent variable are described exactly by 
the regression equation (4). The ratio defined by 


F = S(REG)/(d - 1) (18) 

S(RES)/(n - d) 

follows an F distribution with (d-1) and (n - d) degrees of freedom. This quantity is 
used to determine the statistical significance of the regression equation under consideration 
by comparing it with the appropriate F table value. If the computed F value [equation 
(18)] is greater than the appropriate F table value, the regression equation is statistically 
significant. 


Application 

The application of the preceding analysis to film thickness measurement data is 
concerned with the following equation: 


= b 0 + b 1 X li + b 2 X i li + b 3 X* li (19) 

2 3 

where we wish to determine whether we should use X j j and/or Xjj in the equation. 

This equation can be put in the form: 


Y i " b 0 +b l Z li + b 2 Z 2i + b 3 Z 3i 


where 


Z H = X lr X! 

z 2i = x 2 ir x| 

7 'V' ^ 

Z 3i - Xjj-Xj 
bo = Y 


( 20 ) 


( 21 ) 
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with 



The solution for the regression coefficients in equation (20) (Model 3) then becomes: 


B = S 3 S Y 
3X 1 3X 3 3X1 


where 


B = 
3X 1 


bi 

b 2 

b 3 

Z Z li Z li 
Z Z 2i Z li 

Z Z 3i Z li 


Z Z li Z 2i 
Z Z 2i Z 2i 

Z Z 3i Z 2i 


Z Z li Z 3i 
Z Z 2i Z 3i 
Z Z 3i Z 3i 


S 3 = 
3X3 


(23) 


(24) 


(25) 
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( 26 ) 


xrf 

IZ 2i Y? 

I Y ? _ 

Thus, we can determine the regression coefficients for the following models: 
Model 1: Y® = b^^ Z u 
Model 2: Y® = bg + bj Z u + b 2 Z 2i 

and 

Model 3: Y® = b' Q + bj Z H + b 2 Z 2i + b 3 Z 3i 

These three models yield the following sum of squares due to regression: 

Model 1: S(REG), 

Model 2: S(REG) 2 
and 

Model 3: S(REG) 2 

We now want to test the following hypothesis: 

2 

Hi : X! contributes nothing to variation in Y(i.e., b 2 =0). 

3 

H 2 : Xi contributes nothing to variation in Y(i.e., b 3 =0). 

The test criteria of Hi and H 2 are: 

F ci = , with 1 and n - 3df 

S YY 

F„-> = , with 1 and n - 4df 

C S Y Y 
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where 


S 2>1 - S(REG) 2 -S(REG) 1 
S 3,2 = S(REG) 3 - S(REG) 2 


( 28 ) 


We determine the 95-percent points F[l,n-3] and F[l,n-4] in the F table. If 
F cl <F T , we accept Hj. Otherwise, reject Hi and say that Xi adds significantly to 

the linear regression. If F c2 < Fj, we accept H 2 . If F c2 > F-p, then H 2 is rejected, 
and we can conclude that X? adds significantly to the linear regression. 

RESULTS AMD CONCLUSIONS 


The basic input data for the analysis consist of specimen thickness measurements 
(independent variable) and instrument meter readings (dependent variable). The thickness 
measurements are in mils and the meter readings in counts for the Betascope and voltages 
for the Dermitron. Table 2 is a summary of the six samples from which the basic input 
data were obtained. As seen from this table, the data for a given sample correspond to a 
particular type of coating, base metal, and instrument. A complete discussion of these 
data and the test program from which they were obtained is given in a SPACO internal 
note 2 . 


Linear, quadratic, and cubic regression analyses were performed on the data for 
each sample. As shown in Table 3, the linear models for the six samples were statistically 
significant based on the F test. The hypotheses Hi and H 2 discussed in the section 
Mathematical Development for Multiple Regression Analysis, concerning the quadratic and 
cubic regresson models were then tested. These results are summarized in Table 4. It is 
indicated in this table that for the particular coating, instrument, and base metal used in 
Table 2, the squared and cubic terms postulated in the models did not contribute signifi- 
cantly to the regression. It was, therefore, concluded that the linear regressions in Figures 1 
through 6 are the best choices of models for the indicated data. Confidence limits given 
in these figures establish a range of accuracy at the 95-percent confidence level. 


2. Ibid. 
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TABLE 2. SUMMARY OF INPUT DATA SOURCE 
FOR THE REGRESSION ANALYSIS 


Sample No. 

Instrument 

Coating 

Base Metal 

Number of 
Measurements 

1 

Dermitron 

LB-2 

CB-752 

23 

2 

Dermitron 

SYLCOR 

CB-752 

33 

3 

Dermitron 

SYLCOR 

B-66 

26 

4 

Betascope 

LB-2 

CB-752 

20 

5 

Betascope 

SYLCOR 

CB-752 

33 

6 

Betascope 

SYLCOR 

B-66 

26 


TABLE 3. F VALUES FOR THE LINEAR REGRESSION MODELS 


Sample No. 

n 

F t [ 1, n - 2] , 0.95 
Table Value 

F c , Computed Value 

1 

23 

4.28 

608.4 

' 2 

33 

4.14 

729.9 

3 

26 

4.23 

455.3 

4 

20 

4.35 

87.2 

5 

33 

4.14 

502.4 

6 

26 

4.23 

161.2 
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TABLE 4. F VALUES FOR TESTING HYPOTHESES Hj AND H/ 


Sample No. 

n 

Fj, Table Value 

F c , Computed 

F x [l,n-3] 

F x [l,n-4J 

F cl 

F c2 

1 

23 

4.35 

4.38 

0.02 

0.0001 

2 

33 

4.17 

4.18 

0.001 

0.007 

3 

26 

4.28 

4.30 

0.02 

0.0003 

4 

20 

4.45 

4.49 

0.01 

0.007 

5 

33 

4.17 

4.18 

0.03 

0.00005 

6 

26 

4.28 

4.30 

0.05 

0.002 


a. Hypothesis 
2 

Hi : Xi contributes nothing to variation in Y. 
H 2 : Xi contributes nothing to variation in Y. 
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Y. METER READING 



0 1 2 3 4 5 6 7 

X, INDEPENDENT (coating thickness, mils) 


Figure 1. Linear regression for Dermitron measurements on sample No. 1. 
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Y, METER READING 


X, INDEPENDENT (coating thickness, m X 10‘ 5 ) 



Figure 3. Linear regression for Dermitron measurements on sample No. 3. 
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X, INDEPENDENT (coating thickness, m X 10' 5 ) 







r ♦ - -i ff ■ - f | ■ 

* t { t t M j 


X, INDEPENDENT (coating thickness, mils) 


Figure 4. Linear regression for Betascope measurements on sample No. 4. 
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X, INDEPENDENT (coating thickness, m X 10'®) 


0 2.54 5.08 7.62 10.16 12.70 15.24 17.78 20,32 22.86 



Figure 5. Linear regression for Betascope measurements on sample No. 5. 
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X, INDEPENDENT (coating thickness, m X 10‘ 5 ) 



Figure 6. Linear regression for Betascope measurements on sample No. 6. 
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APPENDIX 


COMPUTER PROGRAM DOCUMENTATION 


This appendix presents operational information on the UNIVAC 1 108 Computer 
Program for the regression analysis application to film thickness measurements experiments. 
The organization of the operational version of the program is depicted in Figure A-l. As 
shov/n, a minimum number of control cards are required to run this deck. A complete 
listing of the program for a typical run is included. 



Figure A-l. Program organization. 
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PROGRAM LISTING 


-RUN ,T REGRES, 3X0390, JUNK I NB I N3 1 3 , 1 .100 

-FOR » I S LQCREG .LQCREG 

C 

C REGRESSION ANALYSIS APPLICATION TO FILM THICKNESS 

C MEASUREMENT EXPERIMENTS 

C 

DIMENSION X ( 100) ,Y( 100) ,YC1( 100 ) , RES 1 ( 1 00 ) , REG1 (100) , TOTAL! 100) »YC 
-2( 100) ,RES2 (100) »REG2( 100) »YC3< 100! ,RES3 ( 100 ) »REG3 ( 100 ) *01(20.2) »D 
-2(20,2 ) »D3< 20*3) »YB1 ( 100) ,YB2 (100) , YB3( 100) , YB4 ( 1 00 ) ., YB5 ( 100 ) ,YB6( 
-100) ,BCDX ( 12 ) , BCDY ( 12) 

DIMENSION DD2<2,2) ,0D3<3,3) ,BJ2(2,2> »BJ3(3,3) 

DIMENSION Z1 (100) ,2? (100) ,23(100) ,YH( 100) 

DIMENSION VARY 1 ( 100) ,VARY2( 100) ,VARY3 ( 1 00 ) ,SGY1 (100) ,SGY2( 100) , 
-SGY3(100) 

DATA (BCDX(I) » I = 1 , 1 2 ) /6HX , I NDE , 6HPENDEN , 6HT ,9*6H / 

DATA (BCDY( I ) ,I = 1»12)/6HY ,11*6H / 

100 FORMAT ( 1 H 1 ) 

101 FORMAT (//) 

102 FORMAT (1H1.74H REGRESSION ANALYSIS APPLICATION TO FILM THICKNESS 
-MEASUREMENT EXPERIMENTS) 

103 FORMAT (64H L I NEAR , QUADRAT I C , AND CUBIC REGRESSION ON INPUT MEASURE 
-MENT DATA) 

104 FORMAT ( 2X , 13HX , I NDEPENDENT , 3X , 1 1HY .DEPENDENT , 3X , 1 2HYC1 , DEGREE-1 , 3 
-X, 12HYC2 > DEGREE-2, 3X.12HYC3, DEGREE-3) 

105 FORMAT •( 2X , F6 . 2 , 1 OX , F8 . 3 ,6X , F8 . 3 , 7X , F8 . 3 , 7X , F8 ,3 ) 

106 FORMAT (31H LINEAR REGRESSION COEFFICIENTS) 

107 FORMAT ( 3X , 3HRO=E13 . 6 , 3X , 3HR 1 =E 13 . 6 , 3X , 8HF VALUE=E13.6) 

108 FORMAT (3X,8HSIGMA Y=E1 3 .6 ,3X ♦ 5HSGR0=E1 3 .6 , 3X , 5HS6R 1 =E1 3 . 6 ) 

109 FORMAT (34H QUADRATIC REGRESSION COEFFICIENTS) 

110 FORMAT (3X,3HSO=E13-«6,3X,3HSl=E13.6,3X,3HS2=E13.6,3X,8HF VALUE=E13 

-. 6 ) 

111 FORMAT ( 3X , 8HS I GMA Y=E 1 3 . 6 , 3X , 5HSGSO=E13 . 6 , 3X , 5HSGS 1=E13 . 6 » 3X , 5HSG 


112 

113 


116 

117 

118 

119 

120 
121 

123 

124 

125 

126 


-S2 = E 1 3 « 6 ) 

FORMAT ( 42H F VALUF FOR SIGNIFICANCE OF X-SQUARE TERM) 

FORMAT (3X,4HF21=E13.6> 

FORMAT ( 30H CUBIC REGRESSION COEFFICIENTS) 

FORMAT (3X,3HTO = E13.6,3X , 3HT 1 =E 1 3 . 6 , 3X , 3HT2 = E 13 .6 , 3 X , 3HT3 = E 13 . 6 » 3X 
- ,8HF VALUE=E13,6) 

FORMAT (3X.8HSIGMA Y=E1 3 .6 , 3X , 5HSGT0=E1 3 ,6 , 3X , 5HSGT 1 =E1 3 . 6 , 3X , 5HSG 
-T2=E13.6»3X,5HSGT3=E13.6) 

FORMAT ( 40H F VALUE FOR SIGNIFICANCE OF X-CUBE TERM) 

FORMAT (3X,4HF32=E13.6) 

FORMAT ( 54H DETERMINANT VALUE FOR QUADRATIC AND CUBIC REGRESSIONS) 
FORMAT (3X,5HDET2=E13.6,3X,5HDET3*E13.6> 

FORMAT (33H INVERSE MATRIX X ORIGINAL MATRIX) 

FORMAT ( 2 1H QUADRATIC REGRESSION) 

FORMAT (17H CUBIC REGRESSION) 

FORMAT (3X»9(E13.6,1X) ) 

FORMAT (13) 

01 (1,1 )=n. 
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PROGRAM LISTING (Continued) 


NAMELIST /INPUT /X s Y * TVAL1 *TVAL2 *TVAL3*NN*XL*XFU YB*YT 
CALL I DENT (335) 

INTEGER P 

READ (5*126) NCASES 
300 WRITE (6*100) 

READ (5* INPUT) 

WRITE ( 6 9 INPUT ) 

INITIALIZATION SUMMATIONS 


BJ = NN 
AX1*0* 

AX2*0a 

AX3*0* 

AY1*0* 

DO 301 I*1#NN 
AX1*X< D+AXl 
AX2*XU)«m+AX2 
AX3=X< I )**3+AX3 
301 AY1~Y ( I ) 4*AY1 

BX1-AX1/BJ 
BX2-AX2/BJ 
BX3*AX3/B J 
BYsAYl/B J 
AVY^BY 

DO 3010 I =1 *NN 

zim«xm-Bxi 

Z2U ) =X ( I )«X(I )-BX2 
Z3U )=X( I ) *#3~BX3 

3010 YHU )=Y( I )-BY 
511=0* 

SI 2 = 0® 

S13=0« 

$21=0* 

522=0® 

$23=0® 

S3J-0* 

$32*0* 

533=0® 

S1Y*0« 

$2Y=0® 

$3Y=0© 

DO 3011 I=1*NN 

sii*zk n*ziu )+sn 

$ 1 2=Z1 ( I ) *Z2 ( I ) +S12 
S13=Z1 ( I )*Z3( I )*S13 
S22=Z2( I )*Z2(I)+-S22 
S23*Z2U)*Z3( I)+S23 
S33=Z3 f I )*Z3 ( I ) +S33 
$1 Y=Z 1 ( I )*YH( I ) +S1 Y 
$2Y=Z2 1 1 ) *YH { I J+S2Y 

3011 S3Y=Z3 ( I ) *YH ( I )*S3Y 
S21*S12 

5 3 1 = S 1 3 
532=523 


21 



PROGRAM LISTING (Continued) 


c 

C LINEAR REGRESSION ANALYSIS 
C 

R 1 = S 1 Y / $ 1 1 
R0=8Y-R1*3X 1 
DO 302 1=1 > NN 

C COMPUTED DEPENDENT VARIABLE (CDV) 

Yd ( I ) =Pn+Rt*x ( 1 ) 

C RE.SIDUAL VARIATIONS (RESV) 

RRS1 ( n = ( Y ( n -YC1 I T 1 1**2 
C REGRESSION VARIATIONS (REGV) 

RFGl (I ) = ( YC 1 ( I ) — AVY ) *#2 
C TOTAL VARIATION (TOTV) 

30 2 TOTAL ( I> = (Y( I ) — A V Y ) * * 2 

SR FS1 =0 ® 

SR^GI =o * 

SY Y = 0 ® 

DO 300 7=1, 

C RESinuAL SUM OF SOS* ( RESSQ ) 

SR ESI sR ESI { I 1+SRFS1 
C REGRESSION SUM OF SOS. { REG SO ) 

SRFGl =REf ; l < T 1 + SRFGl 
C TOTAL SUM OF SOS® (TOT sq } 

3f>3 SYY=TOT AL. ( I l+SYY 

RM 1 =5RFSl / (RJ-2* ) 

RV?=SRF G1 
C F VALUE 

FVALU 1 = RN'2 /RM 1 
C 

C QUADRATIC REGRESSION AMALYSIS 

r 

n?n *i )= c i i 

O'? ( 1 ,?)=S12 
02(2,1 ) = S 1 2 
D2 ( ? , 2 ) =S 2 2 
on 3030 1=1,2 
DO 3030 J=l»2 

3030 002 ( I , J ) =r> 2 ( T , J ) 

CALL GASINV(D2*2tDFT2) 

no 3031 Tsl ,? 

DO 3031 J=I 9? 

3031 BJ2 ( I , J )=•"'» . 
no 3032 M = 1 ,? 

DO 3032 P = 1 * 2 
DO 3032 N -1 *? 

3032 BJ2 ( M *P ) =D2 ( V t\ ) *Do? ( N , P ) +B J2 ( M , P ) 
SI =D2 (1 ,1 ) *S 1 Y+D? ( 1 9 ?)*r-?Y 

S 2 = 02(2,1 ) *S 1 Y +0 2 (■ ? » 2 ) # S 2Y 
SO=R Y—Sl *nx5_s?* r 'X‘> 
on 304 1=1, NM 
C CDV 

YC2 ( I )=S0+S1*X( I )+S?*X( I )*X( I ) 
c RPSV 

RFS2 ( I ) = ( Y ( I ) -YC? ( I ) )**2 
C RFGV 
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PROGRAM LISTING (Continued) 


304 REG2 ( I ) - C YC2 C I >-AVY) *-*2 
SRFS2*0. 

SRFG2»0. 

DO 305 I=1*NN 
C RESSQ 

SR ES2 = RE$2 ( I )+SRES2 
C REGSO 

305 SREG2=REG2 ( I ) +SREG2 
SMI =SRES2 / ( B J-3 « ) 

Sm'2«SREG2/2* 

C F VALUE 

FVALU2*SM2/SM1 

C 

C CUBIC REGRESSION ANALYSIS 
C 

03<1*1)#S11 
D* ( 1 *2 ) =S 1 2 
D3( 1 »3>=S13 
D3<2*1)=S12 
D3 ( 2 * 2 ) -S2.2 
D3 ( 2 * 3 ) =S23 
D3 ( 3 * 1 ) =S 1 3 
D3 ( 3 * 2 ) =S23 
D3 ( 3 *3 ) =$33 
DO 3050 1=1*3 
DO 3050 J=l*3 

3050 DD3 ( I > J )-D3 ( I . J) 

CALL GA$INV(D3939DET3) 

DO 3051 1=1*3 
DO 3051 J = 1 * 3 

3051 BJ3 ( I « J ) =0 ® 

DO 3052 M = 1 * 3 
DO 3052 P = 1 * 3 
DO 3052 N = 1 * 3 

3052 BJ3(M#P)»D3(M*N)*DD3(N*P)+BJ3(M*P) 

T1=D3<1*1 )*S1Y+D3(1 *2)*S2Y+D3(1*3)*S3Y 
T2=D3(2*1 )*S1Y+D3(2*2)*S2Y+D3(293)*S3Y 
T3?D3( 3*1 )*S1Y+D3<3*2 ) *S2Y+D3 ( 3 * 3 ) *S3 Y 
T0=BY-T1*BX1-T2*BX2-T3*BX3 

DO 306 1=1 * N N 

C CDV 

YC3( I >=T0+T1*X(I )+T2*X( I ) *X U )+T3*<X< I )**3) 
C RESV 

RES3 ( I ) = ( Y (I ) -YC3 t I i )**2 
C REGV 

306 REG3 ( I ) = ( YC3 { I )-AVY ) **2 
SRES3 = 0 ® 

SREG3=0. 

DO 307 I =1 * NN 

C RFSSQ 

SRFS3=RES3 ( I ) +SRES3 
C REGSQ 

307 SREG3=REG3 ( I )+SREG3 
TM1=SRES3/ ( BJ-4® ) 

TM?=SREG3/3* 


23 



n n n nnnn nnnn n n n 


PROGRAM LISTING (Continued) 


C F VALUE 

FVALU3=TM2/TM1 

TEST FOR SIGNIFICANCE OF INDIVIDUAL TERMS 

W21=SREG2-SREG1 
W3 2= SR EG 3— SR EG 2 
F?1=W21/SYY 
F32=W32/SYY 

STANDARD DEVIATION OF DEPENDENT VARIABLE FOR LINEAR, 

QUADRATIC, AND CUBIC REGRESSION 

SIGMAl=SQRT(SRESl/(BJ-2. ) ) 

SIGMA2=SORT(SRFS?/(RJ-3.) ) 

SIGMA3=SQRT(SRES3/(BJ-4.) ) 

STANDARD DEVIATIONS FOR LINEAR .QUADRATIC *AND 
CUBIC REGRESSION COEFFICIENTS 

SGRO = S I GM A 1 /SORT ( B J ) 

SGR1=SIGMA1/SQRT(S] 1 ) 

S6S0=SIGMA2/SQRT ( RJ ) 

SGS1 =SI GMA2*SQRT ( D2 ( 1 » 1 ) ) 

SGS2 = S I GMA2*SQRT ( D2 ( 2 , 2 ) ) 

SGTO=S I GM A3 /SORT ( BJ ) 

SGT1=SIGKA3*S0RT(D3(1,1) ) 

SGT2=S IGMA3*SQRT ( D3 t 2 ,2 ) ) 

SGT3=SIGMA3*SQRT ( D3 ( 3 ,3 ) ) 

CONFIDENCE- BANDS FOR L I NEAR .QUADRAT I C AND CUBIC REGRESSIONS 

DO 308 1=1, NN 
ZZ1=Z1 ( I ) **2 
ZZ2-Z2 ( I )**2 
ZZ3=Z3 ( I )**2 

VARY1 ( I )=SGR0*SGR0+SGR1*SGR1*ZZ1 
VARY2 ( I ) =SGS0*SGS0+SGS1*SGS1*ZZ1+S6S2*SGS2*ZZ2 
VARY3 ( I )=SGT0*SGT0+SGT1*SGT1*ZZ1+SGT2*SGT2*ZZ2+SGT3*SGT3*ZZ3 
SGY1 ( I )=SQRT ( VARY] (T> ) 

SGY2 ( I ) =SQRT ( VARY2 ( I ) ) 

SGY3 ( I ) =SQRT ( VARY3 ( I ) ) 

YB1 ( I ) =YC1 ( I ) + TVALl*SGYl ( I ) 

YB2 t I ) =YC1 ( I ) -TVAL1 #SGY1 ( I ) 

YB3 ( I ) =YC2 ( I )+TVAL2*SGY2 ( I ) 

YB4 ( I ) =YC2 ( I ) -T V AL2*S.GY2 ( I ) 

YB5 ( I ) =YC3 ( I )+TVAL3*SGY3( I ) 

308 YB6(I)=YC3( I ) -TV AL3*SGY3 < I ) 

IPAGE=48 
WRITE (6,101) 

WRITE (6,102) 

WRITE (6,101) 

WRITE (6*103) 

WRITE (6,101) 

WRITE (6,104) 
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PROGRAM LISTING (Continued) 


DO 311 I =1 ,NN 

WRITE (6.105) X(I)»Y(I)»YC1(I)»YC2H)»YC3(I) 

I P AGE= I P AGE-1 
IF (IPAGE) 310,309,310 

309 WRITE (6,100) 

WRITE (6,101) 

WRITE (6,104) 

310 CONTINUE 

311 CONTINUE 
WRITE (6,100) 

WRITE (6,106) 

WRITE (6,107) R0,R1,FVALU1 
WRITE (6,108) SIGMAl.SGRO.SGRl 
WRITE (6,101) 

WRITE (6,109) 

WRITE (6,110) SO, SI ,S2»FVALU2 
WRITE (6,111) SIGMA 2 »SGS0*SGS1»SGS2 
WRITF (6,112) 

WRITE (6,113) F21 
WRITE (6,101) 

WRITE (6,114) 

WRITE (6,115) T0,T1 , T2 , T3 , FVALU3 

WRITE (6,116) S IGMA3 , SGT 0,SGT1»SGT2,SGT3 

WRITE (6,117) 

WRITE (6,118) F32 
WRITE (6,119) 

WRITE (6,120) DET2.DET3 
I W=2 
1 = 1 

WRITE (6,101) 

WRITE (6,121) 

WRITE (6,123) 

400 WRITE (6,125) ( BJ2 ( I , J ) » J = 1 ,2 ) 

I W= I W-l 

IF ( IW.EQ.O) GO TO 401 
1 = 1 + 1 
GO TO 400 

401 CONTINUE 
IW=3 

1 = 1 

WRITE (6,101) 

WRITE (6,124) 

402 WRITE' (6,125) ( B J3 ( I , J ) , J=1 , 3 ) 

I W= I W-l 

IF (IW.EO.O) GO TO 403 
1 = 1 + 1 
GO TO 402 

403 CONTINUE 
KK = NN 

CALL QUIK3L ( -1 ,XL , XR , YB » YT ,43 » BCDX , BCDY , KK , X , Y ) 
CALL QU IK3L(0,XL,XR,YB,YT ,35, BCDX. BCDY ,— KK , X , YB1 ) 
CALL QUIK3L < 0 ,XL ,XR , YB . YT » 35 » BCDX.BCDY ,-KK , X , YC1 ) 
CALL QUIK3L ( 0 ,XL ,XR ,YB,YT , 35 ,BCDX , BCDY ,-KK ,X ,YB2 ) 
CALL QUIK3L ( -1 ,XL »XR » YB » YT >43 » BCDX , BCDY ,KK > X . Y ) 
CALL QUIK3L ( 0 ,XL , XR , YB , YT , 35 ,BCDX ,BCDY ,-KK , X » YB3 ) 
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PROGRAM LISTING (Continued) 


CALL QUIK3L{0»XL,XR»YB»YT»35,8CDX»BCDY»-KK»X»YC2) 

CALL QUI K3L (0»XL»XR»YB»YT,35»BCDX»BCDY»-KK»X»YB4) 

CALL QUIK3L < -1 *XL ,XR » YB , YT ,43 > BCDX , BCDY , KK , X , Y ) 

CALL QUIK3L ( 0 ,XL ,XR , YB , YT , 35 »BCDX ,BCDY ,-KK , X * YB5 > 

CALL QUIK3L(0»XL,XR,YB»YT»35,BCDX*BCDY,-K.K»X»YC3) 

CALL QUIK3L ( 0 , XL , XR , YB , YT , 35 »BCDX ,BCDY *-KK , X , YB6 ) 

NCASES=NCASES-1 
IF (NCASES.EQ.O) GO TO 312 
GO TO 300 
312 CONTINUE 

CALL ENDJOB 

STOP 

END 

-XOT 
+ 06 

$ INPUT 

X=. 13*. 32 ,.63,. 76 ,.95 , 1.08, 1.27, 1.90 *2. 2 1,2. 37 ,2.53, 2.66, 2.72. 2. 85 » 
3.16,3.17,3.48,3.79,3.92,4.11 ,4.30,4.48,5.06* 

Y=6. ,19. ,26. ,29. ,30 . ,43 . , 37 . , 57 . , 63 . , 61 . ,70 . , 69. , 70 . ,79. , 78 . , 80. , 83. , 

94. . 90. .95. .100. .100. .100., TVAL1=2.080, TVAL2 = 2.086* TVAL3 = 2.093, NN = 23» 
XL=0.,XR=10.,YB=0.,YT=150.» 

$ 

$ I NPUT 

X=. 38,. 70,. 95, 1.58, 1.70, 1.89, 2. 53, 3. 03,3. 16, 3. 29, 3. 35, 3. 48, 3. 54, 3. 67, 
3.73,3.79,4.05,4.11,4.24,4.36,4.43,4.55,4.62,4.81,5.06,5.19,5.25,5.38, 

5.44.6.01.6.14.6.26.6.96, 

Y=19. ,22. ,23. ,32 . ,31. ,25 . ,40. ,49. ,46. ,53. ,53. ,57. ,54. ,56. ,56. ,53. ,66. , 

58.. 6"1 .,68.,66.,73.,64.,69.,69.,64.,78«,78.,79.,85.,84.,89.,85«, TVAL1 = 2®040, 
TVAL2=2 . 042 , TVAL3=2.045, NN=33, 

XL=0. ,XR=10.,YB=0. »YT=130«» 

$ 

SINPUT 

X=l. 39, 1.58,1 .70,1.89,1.96,2.08,2.21,2.46,2.84,3.03,3.16,3.48,3.60, 

3. 79,4. J.1, 4.43, 4. 46, 4. 62, 4.74,4. 87, 5. .06, 5. 69, 5.82, 6.07, 6. 45,6. 96, 

Y = 29»,28», 23. ,46. ,33. ,36. ,43. ,42. ,50. ,47. ,53. ,59. ,67. ,62. ,64. ,70«, 

74. . 75. .74. .73. .77. .81. .86. ,88. ,93. ,90., TVAL1 = 2.064» TVAL2 = 2.069, TVAL3 = 2.074, 
NN=26 , 

XL = 0. »XR = 10. , YB = 0. ,YT = 15 0. , 

$ 

SINPUT 

X=. 13,, 32,. 63,. 76,. 95, 1.08, 1.27 ,1.90, 2. 2 1,2. 38, 2. 53, 2. 85, 3. 17, 3. 48, 

3.79 ,3. §2,4.11 ,4.30,4.48,5.06, 

Y = 19., 59., 152., 97., 129., 107., 106., 147., 159., 136., 165., 172., 165., 1 72 ,, 

194. . 224. ,237. ,229. ,219. ,218. , TVAL1 = 2.101, TVAL2 = 2.110» TVAL3 = 2.120, NN = 20, 
XL=0. »XR=8« , YB = 0. ,YT = 300. , 

$ 

SINPUT 

X=. 38,. 69,. 95 ,1.58 ,1.70, 1.89,2.53,3.03,3.16,3.29,3.35,3.48,3.54,3.67, 
3.73,3.79,4.05,4.11,4.24,4.36,4.43,4.55,4.62,4.81,5.06,5.19,5.25, 

5.38.5.44.6.01.6.14.6.26.6.96, 

Y=64., 68. ,84., 107., 11 2., 103., 12 4., 130., 131., 134., 150., 142. ,155., 142., 

157.. 149.. 160.. 15 5.. 158.. 169.. 166.. 165.. 162.. 162.. 172.. 166., 179., 172., 

181. . 185. .187. .187. .190., TVAL1=2.040» TVAL2=2.042, TVAL3 = 2.045, NN=33» 

XL=0, ,XR = 10. ,YB = C. ,YT = 240. , 

$ 
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PROGRAM LISTING (Concluded) 


$ I NPIIT 

X±l. 39, 1.58. .1.70,1.89 .1.96 ,2. 08.2. 21, 2.46 ,2.84* 3.03 » 3. 16* 3.48 ,3. 60 » 

3. 79, 4. 11, 4. 43 ,4. 46,4. 62, 4. 74, 4. 87, 5. 06 ,5. 69, 5. 82,6. 07 ,6. 45,6. 96, 

Y = 75.» 71. ,5 0., 116., 94., 97., 10 5., 106., 114., 121., 121. .124.. 139.. 130., 

136. ,148., 149. . 15 5., 151., 160,, 158,, 161., 160., 164., 173., 172., TVAll=2.064, 
TVAL2=2 . 069 , TVAL3=2.074, NN=26, 

XL =0 « , XR= 10, ,YB = 0. , YT = 240. , 

$ 

-PMD.BALE 

-FIN 

-FIN 
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DESCRIPTION OF DATA DECK INPUT PARAMETERS 


The card immediately following the @XQT card is the first card of the input data 
deck. This card specifies the number of cases (NCASES) of data that are to be processed. 
The format is of the form +XX (13) and appears in columns 1 through 3. The information 
between the $ INPUT card and the $ card is associated with a specific set of data and is 
input under the nonexecutable NAMELIST statement. For example, the input statement 
in the program is 

NAMELIST/INPUT/X,Y,TVAL 1 ,TVAL2,TVAL3,NN,XL,XR,YB,YT. 

The forms the input data take on include variable name and subscripted variable. 

In the usage above, X and Y are subscripted arrays, and the remaining variables are simple 
variable names. The specific format of the data can be either integer constants (i.e., +218) 
or real constants (i.e., 1.85921E+00, with or without the E notation). The description of 
the variables in the NAMELIST statement follows. 

X — array containing the values for the independent variable 

Y — array containing the values for the dependent variable 

TVAL1 — value from t-table for NN-(1 + 1) degrees of freedom 

TVAL2 — value from t-table for NN-(2 +1) degrees of freedom 

TVAL3 — value from t-table for NN-(3 + 1) degrees of freedom 

NN — number of data points 

XL — left plot limit for the horizontal X axis 

XR — right plot limit for the horizontal X axis 

YB - bottom plot limit for the vertical Y axis 

YT - top plot limit for the vertical Y axis 
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1108 RUN REQUEST & INSTRUCTIONS 


NAME (LAST & INITIAL) 

Junkin, B. 

OP 

c if 

JOB# PROD » 

310390 

bin# bldg- 

313 4663 

RM# 

A-2 17 

RUN ID 

REGRES . 

RUN 

° F 1- 

EST CPU RUN TIME 
0- HRS, -J .MINS 

CORE SIZE 

PUNCH 

C3 COMPILE 
K) EXECUTE 
a SORT 

Q LANGUAGE max 

ta exec viii F-IV !6*(f 5 

□ SPECIAL FORMS 
TYPE COPIES 


DOES THIS JOB HAVE A RESTART PROCEDURE? □ YES 2D NO 



a OVER 


MtC*ftO FILM 

COPIES 

COPY 

FLO 1 

O P E R ; INIT, 

#FILE$ I 

(■ FRAMES | 

P | 

Fl 

P | 

[ F 

SEQ»# 

1 




: 




OPERATOR COMMENTS: d SEE TECH. □ SEE OPER, 


MSFC - Form 3019 (Rev August 1969) 


□ OVER 


Figure A-2. Job card example. 



TYPICAL OUTPUT RESULTS FOR EXAMPLE NO. 2 


— S'HiPOT- 

X 


Y 


TVALl » 

TVAL2 
TV AIL 3 
NN “ 

XL « 


* 380 QQQQ 0 £+ 00 j 

• X70Q0000E+01 8 
#3l600QQ0£+ai 9 
#35*»QQ000E + Q1 8 

HQSQOQQOE+Ol 9 
4 M 3 ooooqe+oi s 
50400000E+01 s 
SHHOOOOO^ + Ol s 
69600000E+01 s 
OQOOOQOQE+aOs 
OQQOQOOOE+Q 0 9 

ooooooooe+oo* 

QOOQQQOQE+QQs 

oooaboooE+ao » 

QOOQQQQQE + QQ 9 
QOOOOOQOE+OO I 
qoogooqo£+qo S 

OQOOGOdOE+009 

OOOQOOQQE+OQ* 

doU'o'flodoE+obi 

OOQQGOQQE+OO 9 

b-QoobbobE+ooi 

OOOOOOOOE+oo 8 
OOOOGOOOE+OO t 
QOOOQQGOE+OQ 9 

1900000CE+02. 
3 ioooooq£+o 2 9 
H 600 0 0 0 Q E + Q 2 9 
« 5HQ0QQ00E + 02 i 
«T 6 6 00 0 0 0 0 E + 0 2 I 
* 660Q0QGQE + Q 2 » 
® 690G000GE + Q 2 t 
»79000000E+Q2i 
« 8500 Q 00 QE+ 02 i 
•OOQOQOOQE+Od? 
*OOQQaoDQE+ao7 

*OOOOOOOOE+QQ s 
•' OOOOOOOQE+QO 9 
®00000G00E+0D| 
•qooogogqe+oq % 
aOOOQOOQOE + GO * 

®qoooqogqe+qo ? 
®OOOOOOOQE+OQ » 

«OOOOOOQQE+QQ 9 

o qqggooooe+gg s 
*OOQOOGOO£*OO0 
®00000000E+009 
9 OOOOOQGOE+OQ 9 
eOOOOOOOOE + OO 9 

eOOOOOOGGE+'OGa 

e 2OHPOOQOE+0i 

• 2042G000E-5-01 

• 2045QG00E+Q 1 

+ 33 

• OQOGQQGOE + 00 


®7OaO0OQOE*OO 9 
# 18900000E^01 9 
®329OO0OOE*GJ 9 
®367000Q0E+0X 9 
lOOOOOE+Ol 9 
® H5500000E+0 I 9 
*Sl9QQ00GE*0l 9 
96 OIOQOOOE+OI f 
»OGOOQOOOE+O09 
sGOQQOGaQE+QO* 
« OOOOOOOOE+OO 9 
90000GOOOE4*00 9 
•QOOQQOOOE+OOS 
9 OOOOGDGOE + GO 9 
9GQGG00GQE*GQ» 
e ooooqoooe+oo 9 

sOOOODOG.OE + OO 9 
«OOOOOOGGE + O0 9 
tGGQGOOOQE+QQ* 
* QGOQQOQGE + OO « 
«ooQQGaaoE+oo 9 

«0000GD0GE+009 
*00000000 E+OOi 
®QOGQOaQOE + QO s 
®0000QG0GE+aQ8 
*22 QQOGQQE + 02 i 
®2500QQQ0E+02 % 
»53ooqqooe+o2* 
9 56000000E^02 9 
*S8000000E+02 9 
®73QQOGGQE+o2 s 

*6 4 }QG000QE + G2 8 

®850QQQG0E + 02 » 
*OGQOQ 0 QOE+COf 
aQOOOOOOOE+OQ* 
eOOOGGOOGE+OO 9 
900000000E+009 
tQOOCGOOOE+QQ* 
®QCOOQQOCE+OQ 9 
sQOOQGOaOE+OO? 
®ooaaoooQE+oa » 
*0QQ0000QE+00 9 

*GooQoaoaE+aa % 

sQQOOQQOOE+QO 9 
®QCQQQOOQE+CO 9 
sQOOGOOQOE+CO 9 
fflOOOQGOOOE^OO * 
®QG0CG.00Q£*0Q 9 
®QQOOOOQOE+OQ 9 
»GOQOGOOQE+QQ i 


*9500000 OE+OQs 
®25300000E+0l 9 
t 33500Q0QE*Qi 9 
« 37300000E 4>0 1 9 
9 H2H00000E+0 i 9 
®^6200000E^01 9 
*52500000 E+OX 9 
9 6 J ^OQOOOE *0 l 9 
®G000000QE+00 8 
*QOGOOOQQE+Oa a 
9OOOOOOOOE+OO 8 
9 OOOOOOOQE + 00 9 
sQQOQQOOQE+QQs 
9 00000000E+00i 

*ooooooooe+oo 9 

•oaoooaooE+oo* 

®00000000E+009 

«oaooooooE+oo a 

*oooooqooe+oo 9 

* GQQQQQGQE+GQ 9 
*0Q000000E+00| 
®00000000E + 00 8 
eOOOOOOQQE+OO* 
® OOOOOOOQE + 00 9 
*QQQQQQQQE+QQ* 
s230Q0000E+Q2 ? 
®**000000QE+02 j 
8 53000000E+Q2 9 
•560Q0000E+Q2, 

• 61000000E + 02'. 

«5H000000E+02 9 
®780QGQ00E+Q2 , 
®8^000000E+Q2 9 
eOOOOCOOOE+OOs 
sQ000Q00QE+00 9 
9 QQGOOOQQE+GO 9 

,ooooooooe+oo 9 

900000000E+QC* 

®Q00GQ0G0£+G0 9 

®qoooooooe+oo 9 
®ooooqoooe+oq 6 
*QOaaoOQQE+OG 9 
*00000000E + 00 0 
«OOOQOOGOE+OOfl 
*Q0CQQQ0QE+00 0 
®OQQQaOOOE+GQ 9 
*QOQOOOOOE+QO 9 
®Q0QQ00G0E+00 9 
®00000000E+00 9 

*ooooooooe+oo « 


XR 

ss 

. J 0000000 E +02 

— y g- 

* 'SS 

,00000000E + 00 

YT 

& 

. 13000000E+03 

SEND 
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» 1 58Q00QQE + 0 i * 
» 30300000E + 0 i » 
9 3^800000E^01 • 
» 37 9Q0Q0QE *0 1 9 
» H 36Q00QQE + Q i * 
» H8 1 00000E + 0 i * 
e. 53800000 E + 0 X * 
®626Q00G0E+Q1 # 
eOOOOOOOOE+OOs 
eOOOOOOOOE+QO 9 
®OOOOOOOOE+OOi 
®0000000QE+00 9 
®00000000E+Q0 9 

•ooooooooe+oo 9 
*qqgqooqoe+go 8 
*ooooooooE^oa • 

«i QOOOOOQ OE ***00 t 
*OOOOOOOOE+OOt 
oQOOOOOOOE+OO 9 
®QOOOQOOOE+QQ9 
900000000E + 00 9 
9 Q0000000E+00 2 
«OOOOOOQOE+O0 9 
^OOOOOOOOE + OO 9 
oOQOOQOQQE^QQ 
®32000000E + 02 9 

• < *9000000E + Q2* 
9 57000000E + 02 » 

• 53000000E + 02 , 
s^aooooooE^Qz^ 
H A900Q000E+Q2 * 
»78000000E^02 9 
s89000000E+02 8 

®ooooooooe+oo e 
®ooooooooe+oo # 
®oooqooooe*go 9 
^ooooooooe+oo • 

o00000000E+00 8 
sOOOQOOOOE+OO* 
*OOOOOOOOE^OO 8 
sOOOOOOOOE^OOf 
.00000000E+00 9 
9 OOOOOOOOE + OO 9 
sOOOOOOOOE-’-OO® 
«ooooooaoE+oo * 
-OOOOOOOOE+OOt 

«acooooooe+oo « 

® OOOQQOQQE+QQ 9 
«00000000E + 00 9 
« OOOOOOOOE+OO 



TYPICAL OUTPUT RESULTS FOR EXAMPLE NO. 2 (Continued) 


REGRESSION ANALYSIS APPLICATION TO FILM THICKNESS MEASUREMENT EXPERIMENTS 


LINEAR 9 QUADRATJC*AND cubic REGRESSION ON INPUT measurement DATA 


X 9 INDEPENDENT 

• 38 

a 70 

• 95 
1 « 58 
1 *70 
1*89 

_ 2~, 53 ' 

3*03 
3*16 
3*29 
3*35 
3*48 
3*54 
3*67 
3*73 
3*79 
4,05 

h* n 

4*24 
4*36 
4*43 
4*55 
4 ® 6 2 
4*81 
5*06 
5.19 
5*25 
5*38 
5*44 
6*01 
6 » 14 
6*26 
6*96 


y.dependent 

19*000 
22*000 
23*000 
32*000 
31 *000 

25.000 
40*000 
49*000 

46.000 

53.000 

53.000 

57.000 

54.000 
56*000 
56*000 

53.000 

66.000 

58.000 
61 *000 

68.000 
66.000 

73.000 

64.000 

69.000 
69*000 

64.000 

78.000 
'78.000 

79.000 

85.000 

84.000 
.89.000 

85.000 


YC1 $ DEGREE** l 
16*900 
20*646 
23.572 
30.948 
32.352 
34.577 
42.069 
47.922 
49,444 
50.966 
51*668 
53*190 
53.892 
55*414 
56.M7 
56*819 
59.663 
60*565 
62*087 
63*492 
64*311 
65.716 
66*536 
68*760 
71*686 
73*208 
73*91 l 
75*433 
76*135 
82,808 
84,330 
85,734 
93*929 


YC2 # DEGREE-2 
15.516 
19*622 
22 . 8 03 
30.716 
32*207 
34*556 
42.371 
48,370 
49*915 
51,453 
52.161 
53.690 
54.394 
55.914 
56.6X3 
57.311 
60,321 
61,01 1 
62.504 
63.876 
64*674 
66.038 
66*830 
68*974 
71,773 
73*220 
73,885 
75.323 
75.984 
82*200 
83*600 
84*688 
92*291 


YC3 sDEGRE£“3 
19.234 
21*182 
23.058 
28.959 
30.247 
32,377 
40.220 
46*839 
48,600 
50.370 
51*188 
52.963 
53.782 
55*555 
56*371 
57,185 
60*683 
61*481 
63.195 
64,756 
65*656 
67*179 
68.Q54 
70,376 
73,296 
74,745 
75,396 
76,767 
77,380 
82*512 
83,486 
84,313 
87.597 
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TYPICAL OUTPUT RESULTS FOR EXAMPLE NO. 2 (Concluded) 
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